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cn ' Abstract 

, The Sun's l/c^ order gravitational gradient accelerations of Earth and Moon, dependent on 
the motions of the latter bodies, act upon the system's internal angular momentum. Not only 

. does this gravitational spin-orbit force (which plays an important part in calculating gravity 

' wave signal templates from astrophysical sources) slightly accelerate the Earth-Moon system as 

O . a whole, it more robustly perturbs the internal Earth-Moon dynamics with a 5 cm amplitude, 
synodically oscillating range contribution which is presently measured to 4 mm precision by more 

\—{ , than three decades of accumulated lunar laser ranging data. 



1 Introduction 



The Earth-Moon system viewed as a whole is a spinning body moving through the gravitational field of the 
Sun, with the Moon's orbital motion around the Earth contributing five times as much to that system's total 
spin angular momentum as does the Earth's rotation about its axis. In metric theories of gravity, bodies 
with internal angular momentum follow trajectories in external gravitational fields which differ from those of 
so-called test bodies (and simple real bodies); but this difference of trajectories due to a gravitational spin- 
orbit acceleration has not yet been confirmed by observation. The magnitude of the spin-orbit acceleration 
of the Earth- Moon system by the Sun is expected to be about 3^1^ SV/mc^ , and to be radial and repulsive, 
with — GMs/R^ being the strength of the Sun's gravitational gradient, V being the orbital speed of the 
system around the Sun, S and m being the total internal angular momentum and total mass, respectively, 
of the system. The factor 3 in this expression becomes 7-1-2 for scalar-tensor theories of gravity, with 
Eddington parameter 7 being less than one and 1 — 7 being a measure of the fractional participation of 
the scalar interaction in metric gravity. And even this generalized strength parameter is dependent on the 
technical details of defining the spinning body's "center" chosen to express its trajectory (This issue will be 
further discussed later). 
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This spin-orbit acceleration is a consequence of underlying relativistic gravitational gradient forces (accel- 
erations) which the Sun exerts separately on the Earth and Moon; these unequal accelerations thereby affect 
the lunar orbit relative to Earth. The relativistic tidal force acts primarily on the Moon whose internal or- 
bital angular momentum is about 80 times that of the relatively immobile Earth. The lunar orbit's task is to 
transmit the great bulk of this force to the Earth, and it does so by becoming polarized in the solar direction. 
Over 30 years of lunar laser ranging (LLR) has produced data which permits measurement of this orbital 
polarization with realistic precision of 4 mm and with formal error of measurement already below the 1 mm 
level 1^ ; and the continuation of the LLR mission with upgraded technology and improved modeling should 
push this precision of fit for the lunar orbit to even substantially better levels The characteristic size of 
the Moon's acceleration relative to Earth due to the Sun's relativistic gravitational gradient is 'iruVVP' / , 
ru being the Moon's angular momentum per unit mass. This amounts to the fraction 2.6 10~^^ of the Sun's 
Newtonian acceleration, and the frequency of this perturbation on the lunar orbit is synodic. The orbital 
response is the same as would occur due to an Equivalence Principle violation (unequal acceleration of Earth 
and Moon by the Sun) — a polarization of the orbit along the solar direction [Q and with length of 
several centimeters which is an order of magnitude larger than the existing observational precision of that 
polarization's measurement. 



2 The Formal Acceleration of a Spinning Body 

The expected governing tensor equation of motion for a spinning body of mass m and spin tensor J"^) is 
the Papapetrou equation ^ 



whose right hand side will contain further couplings to gradients of the metric field's curvature tensor if 
the body's structure is extended and non-spherical. In Equation (1) m is the body's rest mass, D/dr 
indicates tensor derivative with respect to proper time, = dx^/dr is the body's four-velocity dx^^/dr, 
R'p^y is the metric tensor field's curvature tensor, and with supplementary conditions yet to be specified for 
the spin tensor. This is not a fundamental equation of theory; it must rest on and be derivable from the 
basic gravitational interaction between elements of matter (contributions to a body's total internal angular 
momentum from integer or half-integer Ti quantum mechanical spin of elementary particles are an exception; 
fundamental theory must supply their coupling to gravity). 

A spinning body's motion through a metric gravitational field deviates from the geodesic trajectories 
of ordinary test bodies because of the coupling of the body's internal angular momentum to the curvature 
tensor which is composed of both first and second derivatives of the gravitational metric potentials; therefore 
no fundamental violation of universality of free fall is implied. The leading order acceleration of a body from 
Equation (1), the Papapetrou equation, is 



i, J, fc, ... are spatial indices, o is the time index. (Otherwise in this paper, Greek indices a, /?, ... range 
over all four space and time indices.) Common indices are summed over their range. Assuming that the 
antisymmetric spin tensor J^" is orthogonal to the body's 4-velocity vector J^^V^ = 0, and that its purely 




(1) 
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spatial components are given by the body's internal angular momentum 



gfcy being the completely antisymmetric spatial permutation tensor, e^^^ — 1, then the spin tensor's remain- 
ing components are 

_jok ^ jko ^ y. s)k/c 
In a quasi-static environment of outside bodies, the dominant metric potentials have the form 

goo = 1-2 U/c^ 

U is the outside world's Newtonian gravitational potential, and 7 is one of the PPN Eddington parameters 
which has value one in General Relativity but lesser value in scalar-tensor theories, for example. With 
g = VL/, the curvature tensor components needed above are 

Riook — — [Vg]ifc 

Rtjki ^ -7 {^ki[^g]jk + Sjk[^g]ii - 5ik[Vg]ji - 5ji[Vg]iT^ 

is the 3x3 spatial gravitational gradient tensor of the outside world. The spin- induced acceleration 
can then be written in the form 



\ ((7 + l)[V5]-(l7x5) + 7([Vff]-y)xS') 
a c V / 



with V being the velocity of the spinning system through the outside world's gravitational field. Applying 
this to motion in the Sun's central gravity 

[Vg] = ^^\iRR-\\ 

with ^ = GMg/R^, I being the identity matrix, and R being the unit vector toward the Sun. Neglecting 
the inclination of the lunar orbit plane from the ecliptic, and assuming uniform motion around the Sun, the 
spin-induced acceleration is radial and of strength 

/ SV ^ (m ~ m')m' 
gso = (7 + 2) 2 r!' ^ ru 

with m' being the Moon's mass, while r and u are the Earth-Moon separation and relative speeds, respec- 
tively. As a fraction of the total solar acceleration of the Earth-Moon system, this spin-induced acceleration 
is of size 3 10^^'' and directly unobservable in the foreseeable future. 



3 Relativistic Tidal Force Effects on Motions of Earth and Moon 



Lunar laser ranging (LLR) is here examined as a method for measuring the effect of internal angular mo- 
mentum on the Earth-Moon system's gravitational dynamics; the levels of precision recently reached in the 
on-going LLR mission make this now within reach B. From the fundamental equations of motion for the 
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Moon's coordinate position f{t)i and the Earth's position f{t)2, a "center" coordinate for describing the 
system's dynamics R and another internal coordinate for giving the Moon-Earth separation r can be formed 

mi (l+(v -ui- 55-^1) /c^) n + m2 (l+(v-U2- 55-^2) /c^) r2 

R = (2) 

mi + m2 

f = fi - f2 (3) 

with g being the solar acceleration at the chosen system center R. The inverse of these coordinate definitions 
are also useful 

ri = R + — (1+ (v -112-19 ■X2) /cA f = R + xi (4) 



f2 = R - (^1+ (v -ui- ^5 • ^i) /(^ r = R + X2 (5) 

with m = mi + m2 . These relationships reveal that the total Newtonian gravitational gradient force on the 
system contributes at relativistic order because of the chosen system center 

mixi + m2X2 = K-V-ur (6) 



An internal coordinate p which shows simpler dynamics more closely related to the observable range 
between Earth and Moon is constructed 

r = p ^ [V ■ pv + jg- pp - 'ygp 2) (7) 

mc^ V / 

which generates pertinent differences between these two coordinate accelerations 

= ^ + "^^'c^' • + + iy-'^9)-pv + 2^{V.Vg)-{pv + vp)) (8) 

with dg/dt = V ■ Vg for a static gravitational field; v = dp/dt and a = dv/dt. 



While the internal coordinate p is directly related to the LLR observations, the "center" coordinate is 
primarily useful in an analytic sense and is not observed at any instant, although its longer term evolution 
is a useful first approximation of where the Earth-Moon system is located in the solar system. A sub-system 
of the total system of gravitationally interacting bodies does not have a center of mass-energy because of 
its interaction energy with the total system, so the particular choice for the system's "center" coordinate. 
Equation (2), is motivated by the expression in the Appendix, Equation (xxxxxxx), for the conserved center- 
of-energy coordinate for an entire and isolated system of gravitationally interacting bodies. Because the two 
bodies of the the Earth-Moon system are held together by the relatively elastic gravitational force, itself; 
LLR readily observes the internal perturbations produced by the basic post-Newtonian gravitational gradient 
forces which are responsible for the system's spin-related dynamics. 

If the Sun's acceleration of the Moon from Equation (xxxxxxxx) is expanded to the first tidal order, with 
the Newtonian and 1 /c^ order motional corrections linear in both V and internal velocity Vi collected, there 
results 

{gi)tidai = 0,^ (^{3RR-xi-xi^ 

+ 2(^) —(^j{3RR-p-p)V-v - {l+j){3RR-p-p)-{Vv + vV)^ (9) 
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and the analgous term for the Earth. In the second line of Equation (9), which is aheady order 1/c^, the 
Newtonian approximation xi = pm2/m is used, but xi is retained on the first Une because of the relationship 
given in Equation (6). 



Dynamics of the Earth-Moon System Center 

A combination of the tidal accelerations useful for evaluating the expression in Equation (2) is 

mi{gi)udal + 1^2 {g2)tidal /i„x 

gtidai = [Wj 

m 



mima 2 



(^{2j-l)V -uiSRR-f-r) - (27 + 2)(3M ■ r - r) • {Vu + uV^ (11) 



The equation of motion for the defined "center" can then be assembled and altogether 

+ ■ [W ] • {ur - ru) - 2u ■ [Vg\ ■ rV) (12) 

m c \ / 

Setting u = S X f and averaging over the monthly motion, (ff) = ([I] — SS)/2, this center's equation of 
motion becomes 

^ = 9{R) + -l^{j{VxS). [V5] + (1 + 7)(t? • [W]) X S) (13) 

But it should be noted that the orbital centers of both Earth and Moon do not coincide with the system's 
chosen "center" coordinate R. Prom Equations (4,5) the locations of Earth and Moon are given by 

n = R — V ■ pv -\ p (14) 

mcr \mc/ m 

^ VxS /mi\2 J mi ^ 
r2 = R 2 — V ■ pv p (15) 

Both bodies have a common and slowly varying shift of their orbits' positions in the solar system; and each 
orbit viewed from the solar system rest frame shows the longitudinal oscillation in position due to special 
relativity's time transformation between inertial frames; 

t = V-p/^ + ^l-V^/c'^t' 

which for a given solar system bary-centric time t' , gives the time t for finding the interbody vector at p{t) 
when viewed in the Earth-Moon system's instantaneous rest frame. Defining the actual " center" about which 
the Earth and Moon orbits move 

Rc = R X S 

then requires the adjustments 

(fiR (P Rc r^^i\ ^ 
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in which (fiV/dt^ = V x [Vg] is used, and 

g{R) = g{Ro) + {V x S) ■ [V5] 

The reexpression of Equation (13) in terms of Rc instead of R then shows a spin-orbit acceleration in 
agreement with that which emerges from the Papapetrou equation. 

d'^Rc X , ^ 

= 9{Rc) + 9so 

but with light now shed on what system "center" coordinate achieves that agreement. For a given Earth- 
Moon system orbital frequency f2 about the Sun, the orbital radius is then given by 

3 GMs (, , 1 



(^c)^ = ^ 1 - (7 + 2) 



which indicates a spin-orbit force induced shift in radial location of both Earth and Moon in amount (7 -|- 
2)SVl'imc^ ^ 1.6 mm. 



Dynamics of the Ecirth-Moon separation 



The relativistic solar tidal accelerations of the Moon and Earth, Equation (9), also produces a perturbation 
of the Earth-Moon separation which is directly measured by LLR. 

{g 1) tide - {g2) tide = ^^{iRR-r-f) 

+ "^^rnJ"^ (iV -vimk- p- p) - {-f+\){mR-p-p)-{Vv + vV^ (16) 

which coupled with Equations (8) produces the net internal equation of motion 

d^p Gm _ ^ r^^, 

di = -^p + p-i^9] 

+ ^ ((27 + 2) (V- vp- [Vg] - f • [fg] -vV^ + ■ [Vg] ■ tJf) (17) 
with [V^] = {3RR — I) being the tidal tensor of the Sun's gravity. 

The Newtonian tidal acceleration in Equation (17) produces the classic lunar variation which was first 
discovered analytically by Newton using his universal law of gravity. This perturbation consists of an orbit 
elongation transverse to and a contraction parallel to the solar direction, resulting in a range oscillation with 
half month period. 

f 11 \ 

Spvar - i-cos{'2D)p + —sin{2D)ijj (18) 

with D abeing the Moon's synodic phase (angle between Sun and Moon as seen from Earth). The relativistic 
tidal perturbation given by the second line of Equation (17) perturbs the lunar orbit at the synodic frequency. 
Its projections onto the radial and tangential directions of the lunar orbit are 

gtide -r = (27 + 4) (^^^^^^) ruV cos{D) 

gtide -u = 
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in agreement with a previous analysis. The resulting radial perturbation in the Earth-Moon range is 



c- , . , , f m2 -mi\ f 1 -2n/uj + ...\ Vrir'^ , , , , , , 
&rii) ^ i, + 2) [-^) [ ,_J/^ + J -^cosiD) ^ i6cm)cosiD) 

with the amplifying factor 1/(1 — Til/ui) which almost doubles the perturbation being due to an unusual 
feedback mechanism between this perturbation of synodic frequency and the existing variation disturbance 
of the orbit at twice the synodic frequency [Q. This 5 cm perturbation is well above the present precision, 
4 mm, in measuring the amplitude of the Earth-Moon synodic range oscillation |^[^. This perturbation is, 
of course, part of the outcome which occurs when the equations of motion given by Equation () are computer- 
integrated and the resulting coordinate positions of Earth and Moon are used to determine round trip ranging 
times by light signals. The purpose here has been to make this internal relativistic tidal perturbation explicit 
and show its direct connection with the spin-orbit acceleration of a system as given by Equation (). Without 
this 5 cm amplitude contribution to the Earth-Moon range, the ranging data could not be successfully fit by 
the underlying theoretical model. 
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Appendix 



In scalar-tensor metric theories of gravity, the first post-Newtonian (1/c^ order) gravitational equations 
of motion for bodies i, j, fc = 1, ...N are based upon the lagrangian 
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and the application of a least action principle which produces the dynamical equations 

d dL dL „ ;l TV 

dt dva dra 

The form of the lagrangian given above does not take gravitational self-energy effects within the bodies into 
account; but that lagrangian, or the equations of motion, themselves, can be applied to modify the N-body 
equations for such self-gravity effects. 7 and (3 are two Eddington parameters which are each equal to one in 
pure tensor General Relativity but which generally differ in scalar-tensor gravity . The resulting equations 
of motion are [|| 



2 + 27_ ^, 1 



2/? + 27 ,, ^ \-f- f^H-lrr 1 ^ ii±i2r^ - \ rom 

^2 ~ ( 5^ ( '^J" ~ ^(^J - ^0 ■ + C^u 5j ] (20) 



with the underlying 2-body Newtonian potential and gravitational acceleration of i due to j, and the cumu- 
lative potentials and gravitational accelerations being 

and under the assumption that there are only gravitational accelerations acting among the bodies, 
contributions proportional to accelerations of bodies have been expressed with body accelerations replaced 
by the Newtonian gravity contributions, — > gfe. These equations of motion are computer-integrated 
by different analysis centers to produce the solar system ephemeris for its various bodies. The resulting 
coordinate trajectories supplemented by metric gravity's expressions for the coordinate speed of light 

M ^ c{r^t) ^ (l-(l+7)C/(r-,t)/c2) 
and for the intrinsic rate of generally located and moving clocks 

dr = dt {I- wV2c^ - U{r, t)/c'^) 



8 



are the basic ingredients for modeling the radio and laser ranging data between solar system bodies to Xjc? 
(1/c^ in time) order of accuracy. J7(f, €) is the Newtonian gravitational potential at location r(t) of light 
ray or clock as produced by all solar system sources. 

The equations of motion resulting from this Lagrangian fulfill the classic conservation laws for total system 
momentum, angular momentum, and energy 

dP _ dL _ dE _ ^ 
dt dt dt 

with 

II I 
And the total system possesses a center of energy which does not accelerate 



Applying the complete N-body, first post-Newtonian order gravitational equations of motion given in 
by Equation (), to the Sun-Earth-Moon system, two useful combinations of order acceleration terms 
proportional to the Newtonian accelerations between Earth and Moon, then multiplied by both the Earth- 
Moon system velocity relative to the Sun V and the velocity of the Moon relative to Earth u, are here 
displayed 

(22) 

3 f • f M • f f) (23) 

because in the body of this paper they play an important role in establishing the equations of motion for 
the center and internal coordinates, R and p, respectively. 
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G(m2-mi) 
ai - a2 = 



(2V -uf + V -ruj 

miai+ 171202 Gmim2 , 

= 7 , „ „ [2V ■ ur + 2V ■ ru 

mi + 7712 (mi -|- rn2)c^r'' V 
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